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Abstract 

The competition between spin glass (SG) and antiferromagnetic order (AF) is analyzed 
in two sublattice fermionic Ising models in the presence of a transverse T and a parallel H 
magnetic fields. The exchange interaction follows a Gaussian probability distribution with 
mean —AJo/N and standard deviation Jv/32/iV, but only spins in different sublattices can 
interact. The problem is formulated in a path integral formalism, where the spin operators 
have been expressed as bilinear combinations of Grassmann fields. The results of two fermionic 
models are compared. In the first one, the diagonal S z operator has four states, where two 
eigenvalues vanish (4S model), which are suppressed by a restriction in the two states 2S 
model. The replica symmetry ansatz and the static approximation have been used to obtain 
the free energy. The results are showing in phase diagrams T/J (T is the temperature) versus 
Jo/ J, T/J, and H/ J. When T is increased, Tj (transition temperature to a nonergodic phase) 
reduces and the Neel temperature decreases towards a quantum critical point. The field H 
always destroys AF; however, within a certain range, it favors the frustration. Therefore, 
the presence of both fields, T and H, produces effects that are in competition. The critical 
temperatures are lower for the 4S model and it is less sensitive to the magnetic couplings than 
the 2S model. 

1 Introduction 

There are now several examples of competition between antiferromagnetism (AF) and spin glass 
(SG) in strongly correlated systems as, for instance, heavy fermions and high T c superconductors 
Ej. In some of these systems, there is also a new physics involved with the presence of a 
quantum critical point (QCP) and deviation of the Fermi liquid behavior, so called Non-Fermi 
liquid (NFL) behavior. This raises the necessity of obtaining a solid framework to describe those 
systems where the fermions are operative to give origin to physical processes, such as the Kondo 
effect, in connection with the presence of AF and frustration driven by disorder, particularly, at 
low temperatures where quantum effects become important. 

The complexities of such description can be estimated by the controversies 3 involved in the 
quantum Ising SG which has been investigated by several techniques For instance, one of the 
most interesting open issues is whether the quantum tunneling between the local degenerated min- 
imum of the free energy is able to stabilize the replica symmetric solution of the problem. Recently, 
the fermionic Ising SG in the presence of a transverse magnetic field T has been studied using the 
functional integral formalism [5] . The spins have been represented by bilinear combinations of the 
fcrmion operators. This fermionic problem has been presented in two versions. In the first one, 
the fermionic spin operators have maintained their four natural eigenvalues, where two of them are 
non-magnetic. In the second one, it has been retained only the magnetic ones due to an imposed 
constraint to garantee complete equivalence between the spin and the fermionic problem jj^J. The 
free energy has been obtained in both versions within the replica symmetric theory and the static 
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approximation. It has been found that the freezing temperature Tf(T) decreases to a QCP when 
r enhances in a second order type transition. In that formulation, the replica symmetric solution 
is unstable in the entire SG region. The fermionic formulation is also a natural tool to study the 
interplay between frustration driven by disorder and, for instance, Kondo effect at low temperature 
when quantum tunneling is important. Actually, quite recently, the SG problem with transverse 
field has been studied successfully in the disordered Kondo lattice with this technique [7]. 

The purpose of the present work is to investigate the SG/AF competition in the context of 
the fermionic representation for the Ising spins variables when a magnetic field is applied with 
two components: one parallel (H) and other transverse (T) to the z direction. The component 
r introduces a spin flipping mechanism in the problem which can lead the phase boundaries to a 
QCP [S]. The model used here is the fermionic version of the Korenblit-Shender (KS) model 8 
introduced to study SG/AF competition with classical Ising spins. In the KS model, there are two 
sublattices, but only spins in distinct sublattices are allowed to interact with a random Gaussian 
coupling Jij. The presence of the magnetic field H breaks the symmetry between the sublattices. 
As consequence, it can introduce unusual effects as it has been shown in the mean field theory with 
replica symmetry for the classical interlattice frustrated model 8 . For instance, the field H can 
favor frustration within a certain range. This is in contrast with the well known result from the 
Sherrington-Kirkpatrick model (SK) |9j where the freezing temperature Tf, associated with the 
Almeida- Thouless instability, decreases monotonically with the field H. In fact, the dependence 
of the random internal field with H in the KS model could explain its odd behavior [Hj. The 
internal field applied h p in a particular sublattice (p = a, b) depends on the magnetization m p ' 
and the SG order parameter q p / (p' =^ p). Due to the break of the symmetry between the two 
sublattices, when H is increased in a particular range, one may have a nonmonotonic behavior in 
the order parameters m p i, q p ' and, consequently, in the h p . Therefore, Tt can be enhanced and 
the non-trivial ergodicity breaking region (SG) is enlarged. 

This leads several issues for the SG/AF competition in the fermionic representation of the KS 
model, when both field components (H and T) are applied. Is the scenario described previously for 
the classical KS model preserved or not even for T = case? The answer is not obvious because the 
fermionic representation of the spin operators introduces an important difference as compared with 
its classical counterpart. The replica diagonal spin glass order parameter for the two sublattices 
is not constrained to the unity. It must be solved together with the other order parameters, 
the replica off-diagonal SG order parameter and the magnetization for the two sublattices. As 
consequence, it would introduce a new component in the random internal field. If the previous 
question is answered positively, is there any range of H which can favor frustration when r / 
as in the classical KS? The presence of H, by a mechanism similar to the its classical counterpart, 
could enhance Tf while the T component tends to suppress frustration leading the Tf to a QCP 
0. Therefore, the fields H and T enforce two competing mechanisms in the problem which could 
deeply affect the SG/AF phase boundaries. 

The quantum mechanical partition function of the problem has been obtained following the 
approach introduced in reference Therefore, the functional integral approximation is used to 
deal with non-commutativity of the spin operators which are represented by bilinear combinations 
of Grassmann variables. One important aspect of the fermionic representation of spins Sf is that 
(see Ref. [S]) it has four eingenvalues per site, where two of them are non-magnetic. For the rest 
of the paper, this representation is named 4S model. To recover the usual spin representation, a 
constraint is introduced in order to maintain only the magnetic eigenvalues. This version of the 
problem is named 2S model. The disorder in the problem is treated using the replica trick where 
the order parameters are obtained within the replica symmetric ansatz 9 . It should be noticed 
that the possible occupation of non-magnetic states can produce differences concerning the phase 
boundaries between the two models. In fact, for the one lattice fermionic SG with transverse field 
0, {Tf) 2S > (T/) 4S for T << T c where T c is the value of the field at the QCP. This also arises the 
question how these representations 2S and 4S respond (in terms of the phase boundaries) when H 
and r are turned on. In order to allow us a better comparision between the two representations of 
the problem, the 4S model is kept in the half-filling. 
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In the present approach, the time dependence of the spin-spin correlation functions has not been 
considered (static approximation) |10| . The main argument to justify the use of this approximation 
(see Ref. POE]), even at low temperature, is that our goal is to obtain the phase boundaries in the 
SG/AF competition in the spin fermionic representation and how these boundaries are affected by 
the presence of H and T. As it can be seen elsewhere [HE|, that approximation is quite reasonable 
if the intention is not to reveal the complex nature of phases at very low temperature, but mainly 
to produce phase boundaries which can mimic the experimental ones. 

It is hoped that the results presented here can be a first step to provide a framework which can 
allow us to study the phase boundaries SG/AF competition present in several problems of many 
interacting fermions, such as the heavy fermion system CeiAu\- x Co x Sis |13| where disorder, 
competing RKKY interaction and the Kondo effect have produced a SG alike state and an AF 
phase where the Neel temperature seems to decrease towards a QCP with no trace of NFL behavior. 

This paper has the following structure. In section II, the model is introduced, the saddle point 
free energy and the corresponding order parameters are obtained. In section III, the numerical 
solutions of the order parameters allow us to construct phase diagrams to show the AF/SG com- 
petition. Particularly interesting, for the purposes of the present work, are the phase diagrams in 
the space temperature versus the components of the magnetic fields H and T. It is also shown the 
behavior of the susceptibility. In the last section, the discussion of the results previously shown 
and the final remarks are presented. 



2 Model 

The model considered here is a fermionic Ising model ^1] represented in two sublattices a and b 
where there are two magnetic fields applied: T and H transverse and parallel to the Ising spins, 
respectively. One important point is that only spins located in distinct sublattices can interact like 
in the KS model 0. Thus, 
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where the sums are run over the N sites of each sublattice (a or 6) . The coupling Ji a3h is assumed 
to be a random variable with a Gaussian distribution given by: 



N 

p(J ^ ) = Ve4^ exp 



(^+4J /AQ 2 
64 J 2 



(2) 



The spin operators in Eq. (JJJ are defined as |Sj: 

Sf p = gkt - Ki\ > K = \[c\ v} c lpi + 4 pl a p1 ] (3) 

where hi pCr = c] a c% v a is the number operator, c\ a {c% p a) are fermions creation (destruction) 
operators, with a =] or J, indicating the spin projections, and the sub-index p — a or b represents 
the sublattice. 

In Eq. ©, the spins have been written as bilinear combination of fermion operators which act 
on a space with four states per site (|0 0), |f 0), |0 J.), |T4))- Therefore, Sf has four eigenvalues: 
±1/2 (when there is one fermion in the site i p (a =f or |)) and two when the site i p is unoccupied 
or double occupied (a =f and |). In the present work, two formulations are considered: one 
unrestrained the number of states of the Sf (4S model) , but it considers the average occupation 
of one fermion per site, another restrained the occupation number to n% p -\ + = 1 (2S model). 
This last representation allow us to study the problem avoiding the presence of unoccupied and 
double occupied states |SJ- 

The partition function is given in the Lagrangian path integral formalism where the spin oper- 
ators are represented as anticommuting Grassmann fields (</>, <jf). The partition function for the 
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2S model must consider only states that have one fermion per site. This restriction is obtained by 
using the Kronecker 8 function (S(hi P i +hi P i — 1) = ^ ^ dxi p e M: *»>(™ i P T+ni » > - L_1 )). Therefore, the 
partition function for both models can be represented in a compact form as [5]: 



Z{»} = nilf r dx ip e-<*» f D{(p*4>) exp^}) 
j-. a Jo J 



(4) 



p— a, 6 i p 



where fii p = for the 4S model, which corresponds to the half-filling situation, or fj,i = ixi p for 
the 2S model. The action can be Fourier transformed in time. Thus, we have: 



with 



A{fi = A a M + A\ r + A SG 

= E E £ p M l iuJ + ^ + ^ + Pe?] t p M> 

n i a jb 

s ip (n) = J2<l>t (w + nw M, 



(3 the inverse temperature and the matrices in Eqs. JSHHJ are defined as 
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with the Matsubara's frequencies ui — (2m + l)ir and fl = 2mir (m = 0, ±1, • • • ). In this work, the 
problem is analyzed in the static approximation, which considers only the term when Q = in Eq. 

nanniE]. 

The free energy per site is obtained by using the replica method: f3F = — In Z = — lim n _>o 
1/ 1 (nN)(Z{n) — 1), where Z(n) = (Z n ) j iaj is the configurational averaged replicated partition 
function. The average over P{Ji a j b ) can be performed using the Gaussian distribution given in Eq. 
©.Thus: 



zn = n n n ^ / / e XP e^m 

p=a,b i„ a=l ^ [a=l 
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where a denotes the replica index and S" = S" (0). Eq. (|10|) can be rearranged reviewing the sums 
over different sublattices by square sums over the same sublattice [H] ■ Then, these quadratic terms 
are linearized by using the Hubbard-Stratonovich transformation. This transformation inserts 
auxiliary fields {M^, M£, Q3 7 , and Q^ 7 } on the partition function, therefore: 
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dU exp 
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where / dU = U r =a,b,a R<* dM r JZo U r = a ,b,3 IL 7 d QT' ', and the functional part is expressed 



as: 



e Q7 {//}= IlII^ /" 1 dx^e-^p J D(<f ai 4> a ) 

p—a.b i p ,a ^ 



exp 



E 

p=a,b 



(12) 
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In the thermodynamic limit, the set of integrals in J dU can be performed exactly by the 
steepest descent method, where the auxiliary fields are given by the saddle point solutions: 
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(13) 
(14) 



where (...) denotes the average taken with respect to Eq. 1|12|> . These saddle point equations can 
be used to rewrite Eq. fTTJl as: 



where: 



with 



Z(n) = I ' dU exp | -N [3 2 J 2 £ « 7 " & Jo E «~ E ln ^M 
9£ 7 M = \{^fj dx%e-rt J D[^ p ] exp , 

He / f = E i A M, P ~ Wo»#s?] + 4/3 2 j 2 ]T qp 7s^s;, 



(15) 

(16) 
(17) 



and p = a (p' = b) or p — b (p' = a). 

At this stage, it is assumed the replica symmetric "ansatz" that considers the set g p 7 = q p for 
all a ^ 7, q p a — q p + x P , and m p = m p for all a. The physical quantity \ p = Xp/P where Xp is 
the magnetic susceptibility when Jo = 0. The sums over a in Eq. I|17|) produce quadratic terms 
again, which can be linearized introducing new auxiliary fields in the expression l|16(l . therefore: 



9£ 7 M = / Dz p 

J — GO 

where Dy p = dy p exp (— \/27r, and 



2^ Jo 



K ff 



D[4>* p <j) p ] exp 



The matrix G p 1 {uj) is given by: 

, _ / iuj + n P + /3h p 

Up M = y pr 



0T 

iuj + H P — /3h p 



where 



hp = (H - J m p > + J y/2q p ,z p + J ^J2x p >£, p ) 



(18) 
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with p' p. Therefore, the internal field h p applied in a particular sublattice depends entirely on 
the order parameters of the other sublattice 

The functional integral over the Grassmann variables in Eq. 1)19(1 can be calculated (see Ref. 
[T5]). thus, lnAgj^- = J2uj l n (det G p (w)). In the same equation, the sum over the Matsubara's 
frequencies can be performed like references 114| to give the following expression: 

(22) 



A p e = 2e^ [cosh p, p + cosh (f3y/Ap)] 



where A p = h 2 + T 2 . 

Now, the restriction condition over the number of states can be used for both models. For the 
4S model, the average occupation per site is one. It is obtained putting fj, p = 0, therefore the 
integral over x p is equal to the unity. For the 2S model = ix p ), which corresponds to the spin 
formulation where there are only two states per site (| T,Q), |0,|)), the integral over x p is equal 
to zero. This situation is equivalent to choose an imaginary and temperature-dependent chemical 
potential in the Popov- Fedotov method [T5|. Therefore, the main difference between both models 
is that in the 2S model the contribution of the non-magnetic local states is exactly canceled, while 
in the 4S model it is adopted the half-filling situation. These results are used in Eq. 1(15(1 and the 
free energy can be written as: 



(3F S = (3 2 J 2 (q a xt + qtXa + XaXb) - PJom a m b 



E 

p—a^b 



Dz p In 9 SiP 



In 4 



(23) 



with 
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where s represents the states number allowed in each model. The order parameters q p , \p an( i 
are given by extreme condition of the free energy 1(23(1 : 

JZo D Zp h p sinh^yA^/yA; 



Dz r , 
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e 
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Dz, 



Ho m p [hl cosh(/3VA;)/A p + T 2 sinh(/3 VA;)/(/3A^ /2 )] 



e. 



(25) 
(26) 
(27) 



In particular, when T = 0, the set of integrals over £ p can be performed analytically. In this 
case, it is easy to see that the parameter q p = jj{J2i ^f ^f ) = Qp + Xp — 1 (Eq. (27)) for the 2S 
model. However, in the 4S model, due to the presence of non-magnetic states, q p depends on the 
temperature and the internal field h p . Nevertheless, near T = 0, q p can be expressed by using a 
low-temperature expansion: 



q P = 1 



(s-2) T 
2 JV4o 



• exp 



(g p ) 2 Je~^ 2 / 2 
2 T^/ylf 



(28) 



where g p = 



Jom 



This result suggests that the occupation of non-magnetic states are exponcn- 



J./2q 
V p 

tially small at low T. 

It is well known that the replica symmetric solution can be unstable at low temperature in the 
spin glass phase where the Almeida-Thouless eigenvalue Xat becomes negative. The equation 
for the Xat following as 8 : 



A AT = l-2(/3J) 4 J] f°° Dz p (^- 

p=a,b J -°° ^ S >P 



(29) 



(i 




Figure 1: Phase diagram T/J versus Jo/ J for H = T = 0. The solid lines correspond to the results 
obtained by the restricted model (2S), while the dashed lines correspond to the unrestricted model 
(4S). The dotted lines are the extrapolation carried for lower temperature. Xat = is the Almeida- 
Thouless line. 



where 



f + oo 

i p (z) = e s , p / Dz P 



A 



^ cosh^^/A 



sinh(/3vAp) 



(30) 



3 Results 

The phase diagrams showing the competition between SG and AF for both models with two (2S) 
and four (4S) states can be obtained solving the replica symmetric order parameters given by 
the set of equations (|25l) - (|27|l . In the present work, the numerical solutions for m p , q p , x P an d 
Xat (with p = a or b) are studied by varying the three relevant parameters in the problem, 
which are Jo/ J, T/J and H/J. In this context, the AF phase is given by the order parameters 
l m = m a — TOf, ^ and q a — % =h 0; however, the same situation with the eigenvalue Xat < 
characterizes the mixed phase (SG+AF). It can be considered that, in the SG phase, the order 
parameters l rn = and q a = qb =/= together with Xat < 0. In this quantum fermionic treatment, 
one can identify two situations for the diagonal component of the replica matrix (q p = q P + Xp)- F° r 
the 2S model with T — 0, q p — 1. On the other hand, for the 4S model or for the 2S model when 
r > 0, it is necessary to consider the coupling between diagonal and off-diagonal replica matrix 
elements. Consequently, the parameter x p must be calculated simultaneously with q p and m p , and 
hence it becomes relevant for determination of the remaining physical quantities. In particular, 
for H = 0, the Neel temperature T^o = Tn{Jo,T, H = 0) can be computed expanding the saddle 
point equations in powers of m p . At the second order critical line T/vo, we can make q p = 0, 
therefore X — Xp — ^/{PcJo) and 
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B T ^+r^ 2 cosh/^S+ r2 , 

PcJa{ 7== -Jj - U , 



^ + /z^cosh/W§^- + r2 



(31) 
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Figure 2: Phase diagrams T/J versus T/J for H = 0, (a) < J /J < \[2 (a > 1/V2) and (b) 
Jo/ J = 1.7 (a — 0.59). The line conventions are the same as in Fig. (1). 



where (3 C = 1/T/vo and m = m p = —m p >. For r = 0, the Neel temperature is Tjyo = Jo/[^-^ 
exp ( Two Jo ) ^' wn ^ e f° r -^Vo close to zero, the critical value of the transverse field is given by 
analytical solution of Eq. lpTT|) : Tn = Jq + 2J 2 /J , for both models. 

In Fig. (1), it is shown the diagram T/J (T is the temperature) versus Jo/ J when T and H 
are zero. For that particular situation, the results for 2S and 4S models are qualitative equivalent 
to the KS model with classical Ising spins. Actually, for the 2S model, the Fischer relation 3 
is recovered x a b = X* a bl@ = 1 — 1a,b {x* a b IS ^ ne linear susceptibility for Jo equal to zero). 
Therefore, the 2S model reproduces exactly the KS results |Hj- For high temperature and high 
degree of frustration a = [Jq/J]~ , a paramagnetic phase (PM) is found. The AF solution for 
both models can be found in the region where the degree of frustration a is small enough. In this 
case, m a = —rrib [H]- For increasing a and decreasing temperature, the Almeida-Thouless line in 
Fig. (1) shows the onset (at freezing temperature Tf) of the complex ergodicity breaking region (the 
SG region). There is also an intermediate mixed phase region (SG+AF), where inside the SG, the 
magnetization of the sublattices remains finite with opposite sign. It should be remarked that the 
transition temperatures are different in both models. Particularly, for the same a, (I/)2s > {Tf~)±s, 
{Tn)2S > (Tn)4S and (dT^ /da)2S < (<#Tzv / da)AS {Tn is the Neel temperature). These distinct 
behaviors show the different sensitivity to the magnetic coupling between the two models. For the 
sake of completeness, in Fig. (1), it is also studied the region for Jo/ J < where a ferromagnetic 
(FE) region is found with m a — m&. Thus, one recovers essentially the one lattice spin glass- 
ferromagnetism competition |SJE1- These results show that, in absence of magnetic fields, there 
is a symmetry between the two sublattices. 

If r is turned on (with H = 0), two distict situations can be identified. For the degree of 
frustration a > l/y/2 (0 < Jq/ J < \/2), the magnetization of both sublattices is zero. Therefore, 
the problem is reduced to the one lattice problem in the presence of a transverse field T studied 
in reference |S] (see Fig. (2-a)). For small T, the freezing temperature Tf is different for each 
model. However, when the temperature is decreased and the spin flipping increases due to T, 
the critical behavior of both models tends to become identical. The critical field r c , at the QCP 
for 2S and 4S models, is the same. The Fig. (2-b) shows a phase diagram for a smaller degree of 
frustration a — 0.59 (Jo/ J = 1.7). For this set of parameters, the PM solution is still found at high 
temperature for any value of T. For small values of T, when the temperature is decreased, there 
is a transition to the AF phase at (?V)2S,4S, but (T/v) 2 s > {Tn)4S- For even lower temperature, 
there is other transition to the mixed phase (SG + AF) at (T/) 2 s,4s (again (T/) 2 s > (T/)4s)j and 
finally to the SG phase. These transitions are mainly thermally driven. On the other hand, when 
r is increased, the Neel temperature Tn for both models decreases towards a QCP at T^. The 
temperatures (T/)2S,4S, given by the AT line, decrease when T increases like the Neel temperature. 
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Figure 3: Phase diagrams T/J versus H/J for r = and several values of Jo/ J: (a) Jo/ J = 1.5 
(a = 0.67), (b) Jo/ J = 1.6 (a = 0.63), (c) J /J = 1.7 (a = 0.59), and (d) J /J = 1.9 (a = 0.53). 
The insets represent a zoom in T near Tf. It is used the same convention as Fig. (1) for the line 
types. 



We can observe that both models converge to the same critical behavior when T is increased. The 
reasons for that behavior at a = 0.59 are the same for large degree of frustration (a > l/\/2). 
When we see Eq. I)24[l. at very low temperature, for the 4S model, the term which is integrated 
over £ p is dominant and the distinction between the two models can be neglected. It occurs because 
the contribution of magnetic states dominate the free energy Therefore, there is a distinction 
between the models in the region where the transitions are mostly thermally driven. 

Fig. (3) shows the phase diagrams T / J versus H/ J when r = 0. These phase diagrams are 
studied with four different values of degree of frustration a. The field H breaks the symmetry 
between the sublattices and produces an effect similar to the KS model [HIEj- It tends to destroy 
the AF phase. For instance, the Neel temperatures (Tn)2Sas decrease fast when H is enhanced in 
the region < H < H c (H c is the magnetic field when (Tf) s — (Tjy) s (s = 25, 45)). However, as 
one can see quite clearly in Figs. (3-b)-(3-d), there is a range of magnetic field close to H c where 
the freezing temperatures, associated with the AT instability, increase showing that the frustration 
is favored in this two-sublattice problem [HI . This behavior is different from the one lattice 
problem where the Tf decreases monotonically [§] for any value of H . Nevertheless, for H > H c , 
the behavior of the AT line for 2S and 4S models becomes similar to the one lattice problem. The 
favoring of frustration for both models as function of H can be related to drastic decreasing of 
Neel temperature ^B] near H c . Despite this, the presence of the magnetic field produces different 
effects for 2S and 4S models. The freezing temperature, given by the AT instability, increases 
faster in the 2S model than in the 4S for increasing H. The AF solution is also more robust for 
the 4S model, it is found for greater values of H than the 2S model. 

In Fig. (4), the phase diagrams T/J versus H/J is shown, but with the transverse field tunned 
(T/J — 1.0). The set of values of a is the same as Fig. (3). For that situation, there are very 
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Figure 4: Phase diagrams T/ J versus H / J for V j J = 1.0 and several values of Jo/J: (a) Jo/J = 
(a = 0.67), (b) Jo/J = 1.6 (a = 0.63), (c) J /J = 1.7 (a = 0.59) and (d) J /J = 1.9 (a = 0. 
The same convention as Fig. (3) is used for the line types and the insets. 
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Figure 5: Phase diagrams T/J versus H/J for J /J = 1.7 (a = 0.59) and: (a) T/J = 1.7, (b) 
T/J — 2.4. The same convention as Fig. (1) is used for line types. 



important differences if compared with the T — situation. One effect of the applied T is to 
depress (Tn)2S,4S and (Tf)2S,4S for any value of H. Other important effect is that the AT lines 
for both models have a quite different behavior in the range < H < H c . For instance, in Fig. 
(4-a) where a = 0.67 (Jo/ J — 1.5), when H is enhanced until H c , the Tf, associated with the AT 
instability, decreases faster than the equivalent phase diagram with T = (see Fig. (3-a)). When a 
is increased (see Figs. (4-b)-(4-c)), the AT line shows the delicate balance between the mechanisms 
enforced by T and H. For H « 0, where the presence of T is dominant, the spin flipping tends to 
destroy frustration (dTf/dH < 0). As long as H increases, the frustration is favored again. When 
H k, H c , it is quite clear that dTf/dH > 0. Particularly, close to H c , there is a strong increase 
in the Tf as a decreases. When the T field is enhanced as in Fig. (5) (where T/J — 1.7 and 
T/J — 2.4 with the degree of frustration a = 0.59), the SG character is entirely suppressed. The 
Neel temperature decreases towards H\. The differences between 2S and 4S models still exist for 
H ss 0. However, the increase of T shifts the Neel temperature of both models towards zero. As 
consequence, the difference (T/v)2S — (T^as is increasily small. 

Finally, the study of the magnetic susceptibility allows one to confirm the position of the 
Neel temperatures and to study in detail the differences between the 2S and 4S models near T^. 
The magnetic susceptibility \ is obtained differentiating the saddle point equations (|25l 1261 127|) 
with respect to H. In the limit of zero field (H — > 0), the six- independent saddle point equations 
system is reduced to a three-independent saddle point equations with the following relations: 
m a {H = 0) = -m h {H = 0), q a {H = 0) = q b {H = 0) and x = Xa{H = 0) = X b(H = 0). In this 
case, the resulting linear system for dm p /dH, dq p /dH and dxp/dH (p — a or b) is also simplified 
giving x = dm a /dH = dnib/dH, dq a /dH = —dqb/dH, dxa/dH — —dxb/dH |5]. Thus, 

_ x + n 

X ~T+J (x + Tl) {61) 

where x is given in Eq. (|27|) and II is defined in Eq. (|34|l of the Appendix. 

For high temperature, in the PM phase, the susceptibility has the expression x — x/(T+ JqX)i 
as expected. Now, when the temperature is decreased, x shows a brusque change in its inclination 
at Tjv and a cusp at T g (T g is the temperature below Tf, where m p = m p > = 0), but it does not show 
any anomaly at Tf, as we can see in Fig. Ipjl. For intermediaries values of Jo, the susceptibility 
in the AF phase near Tjv does not decrease with the decrease of temperature (inset of Fig. (Oa)). 
However, it exhibits a slow increase until reaching a maximum value, then it decreases. That is 
because of the large number of frustrated couplings. Nevertheless, strong antiferromagnetic average 
coupling, for example Jo/ J = 2.5, x has a maximum in Tn for both models (Fig. (JJU-c)). In Fig. 
JHJ-b), for Jo/ J = 2.1 and T = H = 0, x exhibits a maximum at Tjv only for the 2S model, while 
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Figure 6: The figures (a), (b) and (c) are plots of \ versus T/J for r = H — and several values of 
J (1.7; 2.1; 2.5). The figure (d) is the x as a function of T/J for H = 0, T — 0.3 J and J = 1.6 J. 
It is used the same convention as Fig. (1) for the line types. The inset in (a) represents a zoom of 
X near the critical points. 



for the 4S model, \ still increases close to T/v in the AF phase. Therefore, due to the presence of 
non-magnetic states in the 4S model, its frustrated couplings are less sensitive than the 2S model 
ones for the increasing of Jo- Another difference of the models is in the value of the critical points 
at high temperatures as showed in Fig. JJjJl. The Fig. (jfjjd) shows the x f° r T = 0.3J, H = 0, 
and a — .625 as a function of T. Again, the critical points are marked by discontinuities in the \ 
that is obtained changing the transverse field. These discontinuities have the same shape described 
above. 



4 Conclusion 

In the present work, the competition between antiferromagnetism (AF) and spin glass (SG) in 
disordered fermionic Ising models has been studied. The Jij coupling among the spins is a random 
gaussian variable with average and variance given by —AJq/N and 32J 2 /iV, respectively. There 
is also a magnetic field applied with components transverse (r) and parallel (H ) to the Ising spin 
direction. The model is similar to the Korenblit-Shender model (KS) [S] for classical Ising spin 
which consists of two sublattices where only the interlattice frustration has been considered. 

The problem has been formulated in the path integral formalism where the spin operators are 
represented by bilinear combinations of Grassmann fields following closely the approach introduced 
in a previous work , which has studied the SG critical properties in presence of a transverse field 
using the static approximation. The focus in this work has been to investigate the boundary 
phase of the SG/AF competition in this two-sublattice model 8 for Ising spins in its fermionic 
representation. For that purpose, the static approximation has been adopted like the previous 
work in Ref. 5^. Two formulations for the fermionic representation for the operator Sf have 
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been assumed: the first one is the unrestricted four-state 4S model (four eingenvalues, two of 
them magnetically insensitive); the second one is the restrict 2S model where the non-magnetic 
eigenvalues are forbidden by an imposed constraint. It should be remarked that the transverse 
component P can tune spin flipping while H destroys the symmetry between the sublattices. 
Therefore, these components have opposite roles as related to the emergence of frustration in the 
problem. 

If r = 0, the Fischer ^Uj relation x p = 1 — g p (p = a, 6) is recovered in the 2S model, as we 
can see from the equations (I25|1 - (I27I) . Thus, the results of the KS model are replicated 8 . That is 
an important difference with the 4S model, where the Fischer relation is not satisfied due to the 
coupling between diagonal and off-diagonal replica matrix elements. Therefore, the internal field 
hp (see Eq. 12 H behaves in a quite distinct way for both models. These differences can be seen very 
clearly in the phase diagram temperature versus Jo/ J given in Fig. (1) with H — 0. The onset 
of AF (at Tjv) and the complex ergodicity breaking given by the Almeida- Thouless (AT) line (at 
Tf) appear for lower temperatures in the 4S model compared with the 2S one. That is a direct 
consequence of the distinct applied internal fields for each model. In the susceptibility x, given in 
Eq. (|32[) . these differences can be seen even in a more detailed way (see figures (JBJ-a)-(JBl-c)). 

When the transverse component T is turned on (H is kept equal to zero) (see Fig. (2)), there is 
an important change in the critical behavior for both models. In the range of degree of frustration 
a = [Jo/ J] 1 > l/\/2 (0 < Jo/ J < V%), the problem is reduced to the one lattice problem studied 
in Ref. [H]- For small T, there are two different transition temperatures from paramagnetism 
(PM) to SG for each model. These transition temperatures converge to same value as T increases 
given origin to the same QCP at T c . For a < l/\/2, for small T, when the temperature is 
lowered, there are the following sequences of transitions: (1) PM/AF at Tjv ((Tn) 2 s ^ {Tn)4s)', (2) 
AF/ (AF+SG), given by the position of the Almeida-Thouless (AT) line at T f ((T/) 2 <? ^ (7»4s), 
where (AF + SG) is a mixed phase; (3) (AF+SG)/SG at T g ((T g ) 2S ^ (T ff ) 4S ), where at this 
transition temperature the sublattice magnetizations become null. For increasing T, the TV for 
each model decreases converging to the same transition line. The same effect is found for Tf. A 
QCP for PM/AF transition is found at T N /J = J /J + 2J/J Q . 

Other important consequence appears when the parallel field H is turned on and, consequently, 
the symmetry between the sublattices is broken. That is source of unusual effects because the 
internal field applied in the sublattice p — a, b depends on the order parameters of the other 
sublattice p' — b, a jSJ^]. For the 2S model, the results of the KS model are reproduced [S] (see 
Fig. (3)). The Neel temperatures for both models (which are different) have a strong decrease as H 
increases. However, the frustration is favored for a particular range of H as shown in Fig. (3) from 
the AT line behavior. It occurs due to the magnetization of a particular sublattice (for instance 
p where m p > m p i ) which contributes to decrease the absolute value of the average internal field 
(hp — H — Jom p i ) that acts on the other one. Hence, in this H range where h p decreases with H , 
Tf can be increased That effect is stronger when a is decreased (Jo/ J is enhanced), since the 
contribution of Jom p > on h p is increased, as shown in Figs. (3-b)-(3-c). Because of the presence of 
non-magnetic state, the 4S model has potentially a weaker response to the magnetic interactions. 
The results show that the AT line can also increase, but slower than the 2S model for the same 
range of H . Nevertheless, in the 4S model, the AF phase is more robust than the 2S one, remaining 
as solution even for higher values of H. However, for very low temperature, the contribution of 
the non-magnetic states in the 4S model becomes less significant. Therefore, both models exhibit 
quite similar critical lines. 

When the transverse field T and H are simultaneously turned on, the effective internal field 



applied in a particular sublattice p becomes A p — ^Jh 2 + T' 2 . Therefore, its effects due to the 

presence of H and T start to compete. The effects from H in the h p are included explicitly, but 
they are also included implicitly by the order parameters. It should be noticed that effects of P in 
the order parameters are present as well. As consequence, there is complex balance of the effects 
discussed in the previous paragraph, dependent on the relation H/T. Particularly, for small H, 
the contribution from T is dominant and the frustration is not favored. As long as H enhances, 
there is a region in the diagram temperature versus H where the frustration becomes favored (see 
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Fig. (4)). Therefore, the transverse component of the magnetic field produces the same effect in 
each sublattice. It suppresses the magnetic order and the frustration. While, in the region where a 
two-sublattice structure is characterized, the parallel component acts in an asymmetric way in each 
sublattice. In this sense, it is always against the AF order. Nevertheless, it can favor frustration 
within a certain range of H . 

To conclude, we studied the SG/AF competition in a two sublattice model [S] where the Ising 
spins have a fermionic representation. There is a magnetic field applied. This field has a transverse 
component (r) which can flip the spins leading the transitions to a QCP. It has also a parallel 
component (H) which breaks the symmetry between the sublattices. These components produce 
opposite effects on the frustration. Our hope is to use this representation to study the several 
strongly correlated problems where there is a SG/AF competition which is associated with a QCP 
as, for instance, in the Cerium alloy Ce 2 Au\- x Co x Si-i |13|. On the other hand, we have used 
the static approximation and replica symmetry ansatz. Quite recently, a new scheme for breaking 
replica symmetry has been proposed j20j which is particularly suitable for the present two-sublattice 
problem studied here. That could be used to improve considerably the description of the SG/AF 
competition of our fermionic model. That would be an object for future work. 
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A Appendix 

In this appendix, the function II (see Eq. is found explicitly: 



where: 



n 

K 2 



/3 2 J 2 K 2 (F 2 -F, 



(3 2 J 2 F 3 



Fi- 

F 2 - 



1 + P 2 J 2 K 1 1 + /? 2 J 2 F 4 

2F 5 (l+[3 2 J 2 K 2 F e ) 



F, 



K 2 {1 + /3 2 J 2 F 6 ) 

l+f3 2 J 2 K! 



l + f3 2 J 2 F 4 
2p 2 J 2 K 2 F 3 F 5 
1 + 13 2 J 2 F 4 



F 2 = 
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F 4 



F, 



Dz 



Dz 



Dz 



dh 2 

d_ 

dh 

d_ 

'dh 



J_ x D^d cosh VA/dh 

s + J^-C^coshVA 

f^DedcaahVE/dh 

s + JZo D t cosh 
IT^D^d 2 coshVA/a/i 2 



Dz 



Dz 



s + f_ x D^ cosh VA 
D£,d 4 coshy/A/dh 4 



I-oo D ^ d2 cosh^A/dh 2 



F 6 = 



oo 32 

Dz- 



s + D£ cosh \/A V s + f-oo D Z cosh ^ 
J^D^d cosh VA/dh d f J^D^d 2 cosh VA/dh 2 * 
s + J^D^coshVA dh \ s + D£ cosh VA 
J^D^d 2 coshVA/9/i 2 " 



dh 2 



with A 



/3 2 r 2 and h = j3h a = /3h b (h a is defined in Eq. J23J) 



(33) 
(34) 
(35) 

(36) 
(37) 
(38) 

(39) 

(40) 
(41) 
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